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Abstract 

Due to the advantage of tunability via size, shape, doping and relatively 
low level of loss and high extent of spatial confinement, graphene quantum 
dots (GQDs) are emerging as an effective way to control light by molecular 
engineering. The collective excitation in GQDs shows both high energy 
plasmon frequency along with frequencies in the terahertz (THz) region 
making these systems powerful materials for photonic technologies. Here, 
we report a systematic study of the linear and nonlinear optical properties 
of large varieties of GQDs ( ~ 400 systems) in size and topology utilizing 
the strengths of both semiempirical and first-principles methods. Our 
detailed study shows how the spectral shift and trends in the optical 
nonlinearity of GQDs depends on their structure, size and shape. Among 
the circular, triangular, stripe, and random shaped GQDs, we find that 
GQDs with inequivalent sublattice atoms always possess lower HOMO-LUMO 
gap, broadband absorption and high nonlinear optical coefficients. Also, 
we hnd that for majority of the GQDs with interesting linear and nonlinear 
optical properties have zigzag edges, although reverse is not always true. 

We strongly believe that our findings can act as guidelines to design GQDs 
in optical parametric oscillators, higher harmonic generators and optical 
modulators. 

Keywords: ZINDO/S, Polyaromatic hydrocarbons (PAHs), Graphene, Clar’s 
rule, Lieb’s theorem. 
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1 Introduction 


Materials with broadband absorption (BBA) and emission, that is, covering 
ultraviolet, visible, and near-infrared regions of the solar spectrum, have important 
applications in photodetectors, broadband modulators, bioimaging, solar cells 
and so forth.^ilf^ Moreover, if the materials with the broadband absorption also 
shows optical nonlinearity they can be very useful in applications involving 
optical parametric oscillation, high harmonic generation,!^ Kerr effectl^ESl and 
multiphoton imaging.^ Thus, finding novel materials with both broadband 
absorption and optical nonlinear activity is of great interest. 

Group IV-VI quantum dots like CdSe, PbSe, CdS, HgTe, ZnSe, etc. have 
already been there in variety of applications involving light emitting diodes, 
bio-imaging, solar cells, and so forth, because of their tunable absorption and 
specific optical nonlinear activity.^i^Hni Materials prepared from high band gap 
semiconductors like ZnS, ZnSe, GaN, and AIN possess ultraviolet optical act ivity 
whereas GdS, rare earth doped GaN materials exhibit near IR activities.! ^^ * ^^ ! 
Although, tuning the size of a quantum dot can vary its active optical range, 
it cannot give the whole range altogether (i.e. simultaneously UV-VIS and IR 
range activity). To this end, GQD s and modified GQDs seems to be promising 
materials for such optical activities.!2212U Together with their higher photostability, 
bio-compatibility and low cost preparation, GQDs may act as a substitute for 
the toxic IV-VI group quantum dots. 

GQDs are the confined graphene materials available in various topologie^211t23l 
and graphene is a layered sp^-bonded carbon material in honeycomb lattice. 
Graphene with its zero band gap has a limitation to its applications in optoelectronics 
due to its zero optical emission. On the other hand, GQDs exhibit a broadband 
absorption and they have emerged as attractive fluorescence materials in the 
ultraviolet, visible and even in infrared regionsDuring recent years, there 
has been a lot of research on the broadband activity of GQDs of different sizes, 
shapes and functionalities through both experiment and theory.!2lH2ll Also, there 
is a progress in identifying the shape and size dependent nonlinear activity of 

GQDs.lsaHaD 

Gonsidering these studies into account, here, we have performed a systematic 
computational study on the linear and nonlinear optical (NLO) properties of 
hydrogen passivated GQDs (hence, may also be termed as polyaromatic hydrocarbons 
(PAHs)) of various sizes, shapes, edge structures and so forth. After careful 
analysis on these GQDs (~ 20) with simultaneous BBA and high NLO coefficients, 
we find that the necessary and sufficient condition for possessing such multifunctionality 
is the presence of inequivalent sublattice atoms. Also, we find that majority of 
the GQDs with only zigzag edges possess this multifunctionality. Additionally, 
we find that some of these GQDs show fascinating 1st hyperpolarizabilities 
10^-10® times larger than the traditional NLO compounds [like p-nitroaniline 
etc]). In the following, first we have described how we have modeled our 
systems and then we have given the details of our computations. Next, we have 
compared the results from our semiempirical calculations on structural stability 
and electronic properties with the earlier studies and then we have presented our 
results on linear and nonlinear optical properties. Finally, we have presented 
the results from first principles calculations on the systems, followed by the 
conclusions and possible extensions to the present work. 
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Figure 1: Some random shaped GQDs with 32 carbon atoms (C32) 


2 Modeling and Computational Details 

As the number of varieties of GQDs which can be generated from graphene are 
huge, following Kuc et. al.l^ we have considered ^ 400 structures, based on 
their size, shape, edge etc. As the hydrogen passivated GQDs have been shown 
to be more stable than the GQDs with bare edges, we have only considered the 
former ones throughout our study. As in ref,l^ we have categorized our GQDs 
as circular (F) or triangular (T) or stripes (i.e. nanoribbons) (S) depending on 
their shape and zigzag (z) or armchair (a) depending on their edges. Thus, F^ 
(T^) represents circular (triangular) GQDs with armchair (zigzag) edges. All 
other GQDs which don’t fit in these categories mainly represent the different 
possible conformers of a GQD with particular number of carbon atoms and we 
refer them as random shaped GQDs. We identify these random shaped GQDs 
with their carbon atom numbers such as G22, G28, C74 etc. In Figure 1 we 
have shown typical examples of random shaped GQDs. 

All the structural optimizations have been performed using self-consistent 
charge (SCC) density functional tight-binding (DFTB) theoryl^ within third 
order expansion of the energy (DFTB3)I^ and with 3ob parameter set,l^ as 
implemented in DFTB-f package.!^ DFTB level of theory is used mainly due 
to the large number of systems (~ 400) considered in this study as well as its 
ability to give trends in band-gaps, energies etc. which are comparable to the 
ones given by DFT, especially for carbon related materials, even with different 
edges, defects and so forth.l^^l^ Geometry optimizations have been performed 
using conjugate gradient method and systems are considered to be optimized 
only when forces on all the atoms are less than 0.0001 Ha/Bohr. As the systems 
are zero-dimensional, we have performed the F-point calculations. For those 
systems whose energy levels near the Fermi-level are almost degenerate, we 
have increased the electronic temperature to 100 K to avoid any convergence 
issues. 

Linear optical properties of all the compounds have been computed at the 
semi-empirical ZINDO/S level of theory as implemented in g09 software package.^ 
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ZINDO/S has been proved to be very successful especially in predicting the 
optical properties of systems containing C, N, O, H atoms like polyaromatic 
hydrocarbon compounds, ^^^^^chlorophylls^^etc.^ At semi-empirical level, nonlinear 
optical (NLO) propert ies of all compounds have been calculated using MOPAC 
2012 program packageAll the first principles calculations for the linear (at 
time dependent density functional theory (TDDFT) level) and nonlinear optical 
properties have been performed using g09. Long range corrected (CAM-B3LYP) 
exchange correlation functional has been used in conjunction with 6-31-|-g(d) 
basis set for all the calculations. A minimum of first 12 lowest excited states 
have been considered in all the studies. GaussSum-2.2.6.ll^ is used to plot the 
absorption spectra and a broadening of 0.333 eV has been used. To ensure the 
reliability of our calculations, we have compared our semi-empirical and first 
principles results on the linear and nonlinear optical properties of p-nitroaniline 
with its reported experimental values. We find that both the results are in close 
agreement (see Table SI of supporting information (SI)). 

3 Results and Discussion 

3.1 Energetic Stability and Electronic Properties: 

All the GQDs considered in this study are found to be thermodynamically 
stable, that is, they have negative formation energy, Eporm = ^tot - 
- Nc*Ec, where Ejoj, Ep and Ec are the total energy of the system, energy 
of the hydrogen atom in a H 2 molecule (i.e. E^f 2 / 2 ) and energy of the carbon 
atom in a graphene lattice (i.e. Ecraph/^c), respectively. Here, Nc and Np 
are the number of carbon and hydrogen atoms in the system. At DFTB3 level 
of theory, we find Ep and Ec to be -9.123 and -44.291 eV, respectively. A 
plot of formation energy per atom vs N///(N//-|-Nc) of all the systems is given 
in Fig. 2a. Clearly, there is a near linear relationship between the formation 
energy per atom and the number of edge atoms in all the systems (notice the 
linear fit in Fig. 2a), that is, system with lesser number of edge atoms is easier 
to form and vice-versa, as expected.^^^^^^^^ Simil ar res ults have been observed 
in some of the earlier studies on GQDs and PAH^^UMI jn agreement with these 
previous studies, we also find that among the different GQD shapes studied 
here, circular GQDs are the most stable ones and ribbon like GQDs are the 
least stable. All other GQDs’ (triangular, random etc.) stability fall in between 
these two types of GQDs [see Fig. 2a]. The reason for such a trend is again 
due to the less number of edge atoms in circular GQDs than in other GQDs 
considered in this study, as evident from the x-axis of Fig. 2a. Recent molecular 
dynamics simulations have also shown that among the different GQDs, circular 
and triangular GQDs with zigzag edges as the most stable ones till ~ 4000 K.l^ 
Next, the energies of HOMO, LUMO and their difference (i.e. HOMO-LUMO 
gap (HLG)) of all the GQDs are plotted in Fig. 2b as a function of number of 
carbon atoms. The calculated HLG values are mainly in the range of ~ 0-3 eV. 

Also, from Fig. 2b and 2c, it can be observed that for a particular Nc, one can 
tune the HLG from ^ 0-3 eV depending on the shape and edges of the GQD. 

Also it should be noted that such tuning is possible even for the systems with 
Nc between 20 to 50. In fact, synthesis of GQDs (actually, PAHs) of different 
sizes have already been carried out.^ From Fig. 2c, it can be noticed that 
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Figure 2: (a) A plot oiT!^Form per atom versus number of edge atoms to the total 
number of atoms (N) of all the GQDs. Straight line shows the linear fit. (b) 
Energies of HOMO, LUMO and the HOMO-LUMO gap (HLG) of all the GQDs 
and (c) Ghanges in the HLG with size for different shaped GQDs. Symbols 
T, F, S represents triangular, circular and striped GQDs. Subscripts a and z 
represents armchair and zigzag edges. See the “Modeling” for further details. 


HLG of the systems with zigzag edges converge rapidly to zero (reaching the 
semi-metallicity of graphene) than the armchair ones, irrespective of the shapes 
and the calculated trend of convergence is T^-GQDs > Sj-GQDs > F^-GQDs 
> Ta-GQDs ^ Fq-GQDs > Sa-GQDs. As HLG reflects the kinetic stability 
of a system, the above trends suggest that kinetic stability will be highest for 
Sa-GQDs and least for T^-GQDs and S^-GQDs. As suggested by the Glars 
rule,l^ higher kinetic stability of Sa-GQDs, compared to the other structures 
is due to the presence of larger number resonant sextets in these structures. 
Similar reasons are also known for the lesser stability of zigzag edged structures 
compared to the armchair ones.^^ One may also notice that the HLG of “Sz 
and Tz”, “Ta and Fa”-GQDs follows similar trend as Nc increases (for Nc > 
60) as has also been observed in some of the recent studies.!^ Finally, as the 
HLG of these GQDs are tunable over a wide range and as HLG can be used as 
a rough estimate for the optical gap,l^ one may immediately expect that the 
optical properties of these GQDs can also be tuned over a wide range and the 
results of the respective calculations are given below. 

3.2 Optical Properties: 

First, we present the optical absorption of all the systems calculated at the 
ZINDO/S level of theory. Here, we have analyzed only the 20 low energy 
singlet excitations from the ZINDO/S results. Absorption spectra of PAHs 
mainly consists of 3 bands, namely, alpha (a), beta (/?) and para (p), out of 
which the most intense ones being /? and p-bands (notations are according to 
Glar’s rule,l^ where p(,d)-bands corresponds to the bands at higher (lower) 
wavelengths). Interestingly, in a very recent study®! ji; jj^s been concluded 
that ZINDO/S is good at predicting the most intense p and /3 bands of all 
G 32 H 16 benzenoid PAHs. Considering these facts, first we have plotted the 
histograms of “wavelengths corresponding to the most intense p-bands (Pmax) 
and /3-bands {Pmax)”, respectively, in Figs. 3a and 3b and the corresponding 
oscillator strengths (OS) histograms in Figs. Sla and Sib. From these figures 
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Figure 3: Histograms of wavelengths corresponding to (a) /3max, (b) Pmax 
excitation in all GQDs. Pmax excitations above 2500 nm have been omitted 
for clarity. 


it can be noticed that, majority of the systems have their /3max and p^ax ia 
the UV-VIS region (200-760 nm) and the oscillator strength of Pmax (Pmax) is 
almost always (for majority of structures) > 0.5 (0.1). Thus, majority of GQDs 
considered in this study absorb strongly in the UV-VIS region (in particular, 
their Pmax (Pmax) i® located in the region between 250-450 (300-700) nm). 

However, interestingly, we find ~ 70 GQDs whose Pmax i® ia IR-region (> 760 
nm). Materials absorbing in IR region are of great interest in the preparation 
of solar cells because half of the solar energy received by earth is in IR radiation 
range and most of the present day solar cells do not utilize this energy region.l^ 
Thus, knowing the reason for the IR-activity of these GQDs will be of great 
use and for this we have analyzed their Pmax transition. In Table SI, we 
have given the calculated wavelength, OS and the major contributions of the 
molecular orbitals corresponding to the Pmax transition for all the IR-active 
GQDs. Clearly, Pmax transition always has the major contributions from excitations 
involving the frontier orbitals (that is, HOMO-I, HOMO, LUMO and LUMO-I-I), 
especially from HOMO and LUMO. Thus, the changes in these frontier MOs 
lead to changes in the Pmax transition. Also, some of the earlier studies on 
PAHs have found that HLG of these systems is almost equal to the energy 
corresponding to the Pmax transition (see ReP^ and references there in). 

For a few of these GQDs, we find HLG to be very small. In general, small 
HLGs occur either due to extended delocalization (as in conjugated carbon 
cha ins) or i f there exists lesser number of resonant sextets (according to Clar’s 
j^uiglUMIinij In our case, however, the very small HLGs are seen due to completely 
different reasons. If we look at the structures of these GQDs closely, we find 
that, they don’t have same number of sublattice atoms (i.e. ^ 0). 

In fact, in all the random shaped GQDs, we find there exists two additional 
sublattice atoms of one type (i.e. | - Nb | = 2). p^ orbitals of these 

additional atoms remains as non-bonding orbitals and appear at the zero of 
energy (i.e. at the Fermi-level) in the energy level diagram. If there were no 
interactions (as in tight-binding calculations), both of these levels would be 
degenerate and would appear exactly at the zero of energy, (similar to what has 
been observed in triangular GQD^^. However, because of interaction terms in 
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Figure 4: Schematic diagrams of four structural isomers of C 32 H 18 GQD are 
given in (a)-(d) and their HOMO isosurfaces are given in (e)-(h), respectively. 
Iso-value of 0.02 e/A^ is used for all the plots. 


ZINDO/S Hamiltonian, we find the two levels to appear above and below the 
zero of energy with a very low energy gap (few meV). Interestingly, these two 
levels have opposite parity, due to which the transition dipole moment between 
the two become non-zero. Thus, these two levels give rise to optical transition 
with a finite oscillator strength (OS). Since, the energy gap between these two 
levels is too small, the optical absorption appears in IR-region. 

As an example, in Figs. 4a-4d we have given four structural isomers (here 
after, addressed as 4a, 4b, 4c and 4d, respectively) of C 32 H 18 , where only 4b 
and 4d have the sublattice imbalance. As explained, only for 4b and 4d, we find 
Pmax ill IR-region (> 2000 nm) but not for 4a and 4c. In Fig. 4, we have also 
given the conjugation and isosurface plots of HOMO for these GQDs. As can be 
seen, because of sublattice imbalance, the conjugation in 4b and 4d GQDs is not 
continuous and there are “conjugation breaks”, which are clear demonstration 
of solitonic structure (conjugated system), domain walls (seen in ferromagnetic 
metal blocks). The main point is that, these defect states are intrinsic in these 
GQDs and these have not been externally induced. 

IR-activity of GQDs which absorb below 2500 nm is mainly due to the zigzag 
edge nature of these GQDs, which lowers their HLG. For example, it is well 
known that the polyacenes have the lowest HLG among the various PAHs^^ 
and Pmax of hexacene (6 fused benzne rings) itself is 750 nm. Also, through 
TDDFT calculations, previously our group has shown that Pj^ax of rectangular 
GQDs is ~ 1900 nm.l^ Inspecting the structures of GQDs which absorb in 
the region of 760-2500 nm, we find that all these GQDs have either polyacene 
type structure or rectangular type structure, with some of their edges being 
substituted with ethene, propene, cis-1,3-dibutene etc. (see Fig. S3) Also, it 
is important to mention that HOMO of all the IR-active structures is different 
from that of the non-IR-active structures. 

Frontier MOs of IR-active GQDs have larger number of nodes, and hence, 
look like the collection of p^ orbitals on individual carbon atoms without overlap 
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(for example, see Figs. 4f and 4h). The reverse is true for the non-IR-active 
GQDs (see Figs. (4e, 4g)). Presence of large number of nodes destabilizes 
HOMO compared to its structural isomers with less number of nodes, and hence, 
lesser HLG and IR-activity. As an example of the above mentioned observations, 
we have given absorption spectra and isosurfaces of HOMO of G32 and C74 
GQDs in Figs. S2-S4. Finally, as the OS of p^ax peak for majority of these 
GQDs is > 0.5 and as OS of /3max peak is almost always found to > 0.5, we 
find that these GQDs can have broad band absorption (BBA), as predicted 
earlier for the rectangular GQDs.^^ BBA of these GQDs can also be seen in 
Figs. S2 and S4. To conclude all the above results, we find that GQDs with 
inequivalent sublattice atoms or GQDs with rectangular or stripe shapes can 
absorb in IR-region and they may be suitable candidates for BBA. 

To put the results obtained from the ZINDO/S method in a solid footing, 
we have performed TDDFT calculations at CAMB3LYP/6-31+g(d) level of 
theory on a few GQDs. First, we will present our results on GQDs of various 
shapes. In Figs. 5a-5f, we have given the absorption spectra of Sq, S^, F^, Fj,, 

Tq and T^, respectively, calculated at both ZINDO/S and TDDFT levels of 
theory along with the iso-surfaces of their HOMO (only from TDDFT). Clearly, 
absorption profiles of both the methods compares farily well, although OS values 
predicted by ZINDO/S are higher than that of TDDFT. Also, Xmax predicted by 
ZINDO/S is consistently red-shifted compared to the TDDFT predicted values. 
Consistent with the previous arguments on the isosurface of HOMO (calculated 
using ZINDO/S), even with TDDFT we find larger number of nodes (see Fig. 

5f) in the HOMO if the GQD has IR-activity and it has more overlapping 
character if the GQD is not IR-active (see Fig. 5a-5e). Also, we find that 
the character of the p^ax excitation (i.e. MOs involved in the excitation) are 
similar in both the methods. Importantly, we find that GQDs whose HOMO is 
mainly localized on the edge atoms (as in S^ and T^) and whose p^ax excitation 
has major contribution from HOMO to LUMO, are IR-active. From Fig. 5, 
one may also infer that the presence of zigzag edges is only a necessary, but 
not a sufficient condition (example being the F^ GQDs) for the IR-absorption. 
Finally, to see the effect of inequivalent sublattice atoms on the IR-activity, we 
have considered five C28-GQDs with N^ - Ns = 2, and we find all of them to 
be IR active, again consistent with the ZINDO/S results (see Table S2). Thus, 
we find that, results of ZINDO/S and TDDFT are consistent and compare well 
for the GQDs considered in this study. 

3.3 Nonlinear Optical Properties: 

In this subsection, we present the linear polarizability (a) and first hyperpolarizability 
(/3) of all the GQDs calculated using the finite field approach as implemented 
in the MOPAC and g09 packages. Expressions for the dipole moment and the 
energy of a molecule interacting with an external electric field are given by Eqns. 

1 and 2.1^ 

Mz — “f ^ij^j 4” 4“ 4- . . . (1) 

~ ~ ~ '^'YijkiFjFkFi — ... (2) 
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E{F) = E(0) - n,F, 



Wavelength (nm) Wavelength (nm) 


Figure 5: Absorption profiles of GQDs of various shapes calculated at both 
ZINDO/S level of theory and using TDDFT at CAMB3LYP/6-31+g(d) level 
of theory. Insets in each figure show the isosurface of the HOMO of that GQD 
calculated using TDDFT. (a)-(f) represents S^, S^, Fq, F^, Tq and T^ GQDs, 
respectively. Iso-value of 0.02 e/A^ is used for all the plots. 
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Figure 6: Histograms of average values of (a) polarizability, a and (b) first 
hyperpolarizability, (3 of all GQDs. 


where, /tq is the permanent dipole moment, ajj, Pijk and 'jijki are the linear 
polarizability, 1®* and 2"'^ hyperpolarizability tensor elements, respectively. Also, 
for a molecule, the avergae values of above quantities {pav etc.) are defined as 


where. 


A 


3 

5 


l^av — 

1 

^av — 2 

i^av — 


(Aii + P 


(Mx + My + 

(3) 

“1“ ^yy O^zz) 

(4) 

(/? 2 +/ 32 +/ 32 ) 1/2 

(5) 

ijj +/3ikk),hj,k = x,y,z 

(6) 


In Figs. 6a and 6b, we have plotted the distribution of average a and /3 
values for all the GQDs at static field. Similar to absorption profiles, majority 
of the GQDs’ a and (3 values are confined to a small region. For these majority 
GQDs, we find that the a and (3 values are in the range of 250—700 a.u. 
40—100 and 1—200 a.u. (10“^^ —10“^° esu), respectively. Gompared to the 
a and /? values of para-nitroaniline (16.346 a.u. and 978.21 a.u., respectively), 
it is nice to notice that majority of the GQDs already have high polarizability 
and reasonable hyperpolarizabilities. Importantly, we find that several GQDs 
possess a and (3 values which are orders of magnitude greater than that of 
para-nitroaniline (see Fig 6a and inset of Fig 6b). 

In general, both linear polarizability and first order hyperpolarizabilities 
have an inverse relationship with the energy gap between the states involved in 
the polarization, and are directly proportional to the transition moment. Thus, 
we can expect an increase in a and j3 if the ground and excited states are closely 
spaced or the transition moment between the states is high or both. From the 
above reasoning, one can also infer that GQDs with low HLG and whose 1®* 
excited state has major contribution from HOMO to LUMO transition should 
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give higher a and /3 values. Indeed, we find that all the GQDs which are 
IR-active also have high a and /3 (except the GQDs with inversion symmetry) 
values, that is, above the range of 250—700 a.u. and 1—200 a.u., respectively. 
Also, we find that some of the GQDs with zigzag edges, like F^, which are 
not IR-active but have very high oscillator strength for the (see Fig. 5d) 
excitation also show higher a values. However, to the presence of inversion 
symmetry, such GQDs do not have higher f3 values. Finally, we again find 
that trends in our results from semiempirical calculations compares farily well 
with that of first-principle calculations (see Table S4). Based on all the above 
results, we conjecture that GQDs with very low HLGs can have both broad band 
absorption and nonlinear optical activity, and hence, are potential candidates 
for optoelectronic devices. 

4 Conclusions 

We have performed a systematic study on the GQDs of various sizes, shapes 
and edges to explore their linear and nonlinear optical properties. First, we 
find the formation energies of GQDs have a near linear dependence on their 
number of edge atoms and HOMO-LUMO gaps of a GQD with a particular 
number of carbon atoms can be tuned from ~ 0-3 eV depending on its shape 
and edge nature. Trends in the HLG can be understood based on the Glar’s rule 
of aromatic sextets for majority of the sytems. Extremely low HLGs of certain 
GQDs is due to the presence of unequal number of sublattice atoms in these 
GQDs, that is, ^ 0. Tunability of HLG has also been reflected in the 

tunability of the absorption profiles in these GQDs. We find that majority of 
the GQDs absorb strongly in the UV-VIS region with their Pmax (Pmax) being 
located in the region between 250-450 (300-700) nm) and their a and (3 values 
are in the range of 250—700 a.u. and 1—200 a.u., respectively. However, ~ 70 
GQDs have their p^a.x IR-region and have higher a (> 700 a.u.) and j3 (> 
200 a.u.) values. A common feature which we find in all these IR-active GQDs 
is the existence of larger number of nodes in the isosurface of HOMO which 
leads to an increment in HOMO energy, and hence, decrement in the HLG. Due 
to their high oscillator strengths in both UV-VIS and IR-regions these GQDs 
can possess broad band absorption. With their high a and /? values along with 
the BBA, we expect them to be potential candidates for optoelectronic devices. 
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(a) Oscillator Strength 
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Figure SI: Histograms of oscillator strengths corresponding to (a) Pmax, (b) 
Pmax excitation in all GQDs. 
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Table Sl:Reported values of the absorption, polarizability and 1st 
hyperpolarizability (at static field) for p-niitro aniline along with the results 
from the present work 


Methods 

Absorption (nm) 

aav (a.u.) 

0av (a.u.) 

Experimental 

350 (in dioxane)“ 

375 (in water) 

1147® 

1072(+/-44)® 

DFT reported 

29T^ 

102.8® 

1794.7® (/?^) 

QM/MM 

2667 

110.33 

978.23 

DFTK 

269 

94.9 

3741 

SemiempiricaD 

316 

100.5 

492.5 


a) Ultrafast Charge-Transfer Dynamics: Studies of p-Nitroaniline in Water and Dioxane, C. L. Thomsen, J. 
Thgersen, and S. R. Keiding J. Phys. Chem, A 1998, 102, 1062-1067. 

b) Experimental Investigations of Organic Molecular Nonlinear Optical Polarizabilities. 1. Methods and Results 
on Benzene and Stilbene Derivatives Ching, L-T.; Tam, W.; Stevenson, S. H.; Meredith, G.; Rikken, G.; Marder, 

S. R. J Phys Chem 1991, 95, 10631. 

c) A comparison of molecular hyperpolarizabilities from gas and liquid phase measurements Kaatz, P.; Donley, E. 
A.; Shelton, D. P. J Chem Phys 1998,108, 849. 

d) Relaxation of Optically Excited p-Nitroaniline: Semiempirical Quantum-Chemical Calculations Compared to 
Femtosecond Experimental Results, Vadim M. Farztdinov, Roland Schanz, Sergey A. Kovalenko, and Nikolaus P. 
Ernsting J. Phys. Chem. A 2000, 104, 11486-11496. 

e) B3LYP Study of the Dipole Moment and the Static Dipole (Hyper)Polarizabilities of para-Nitroaniline in Gas 
Phase. Int. J, Quan, Chem. 2006, 106, 11301137. 

f) Solvent Effects on the Electronic Transitions of p-Nitroaniline: A QM/EFP Study, Dmytro Kosenkov and 
Lyudmila V. Slipchenko, J. Phys. Chem. A 2011, 115, 392401. 

g) The first hyperpolarizability of p-nitroaniline in 1,4-dioxane; A quantum mechanical/molecular mechanics study, 
Lasse Jensen, Piet Th. van Duijnen, J. Chem. Phys. 2005, 123, 074307. 

h) This work. DFT at CAMB3LYP/6-31-|-g(d) level of theory. Semiempirical results of absroption are calculated 
with ZINDO/S and a and 0 with MOPAC. 
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Systems 

Wavelength (nm) 

Osc. Str. 

MO contributions 

c38hf 

764.6684 

0.9090 

HOMO^LUMO (93%) 

c42hd 

781.5389 

0.2922 

HOMO^LUMO (93%) 

c74ha 

784.6053 

9.4406 

HOMO^L+1 (45%), H-l^LUMO (46%) 

c60hm 

799.0677 

1.6014 

HOMO^LUMO (91%) 

c54he 

803.7818 

0.8854 

HOMO^LUMO (91%) 

c78h42 

805.0343 

0.5534 

H-l^L+1 (69%), H-2^L+2 (22%) 

c38ha 

808.7628 

0.3008 

HOMO^LUMO (93%) 

c60hk 

812.7390 

1.1520 

HOMO^LUMO (93%) 

c54hh 

823.9738 

0.6357 

HOMO^LUMO (93%) 

c82h44 

827.3278 

0.5781 

H-l^L+1 (66%), H-2^L+2 (25%) 

c42hh 

827.8249 

0.9441 

HOMO^LUMO (93%) 

c50hd 

856.1794 

0.9842 

HOMO^LUMO (93%) 

c54hc 

859.8012 

0.3036 

HOMO^LUMO (91%) 

c42hc 

860.1591 

0.3094 

HOMO^LUMO (92%) 

c50hb 

882.9464 

0.3115 

HOMO^LUMO (91%) 

c60hb 

893.8313 

0.7016 

HOMO^LUMO (82%) 

c54hg 

900.3873 

0.9397 

HOMO^LUMO (93%) 

c60ho 

902.2876 

1.7710 

HOMO^LUMO (94%) 

c60hf 

904.4597 

0.6102 

HOMO^LUMO (92%) 

c74hs 

905.2522 

0.0556 

HOMO^L+1 (72%), HOMO^L+5 (11%) 

c60hh 

908.3694 

1.0712 

HOMO^LUMO (92%) 

c60hp 

910.4372 

1.8297 

HOMO^LUMO (94%) 

c60hr 

920.8507 

1.8023 

HOMO^LUMO (94%) 

c54hb 

923.8699 

0.3222 

HOMO^LUMO (90%) 

c74hr 

928.7141 

1.1620 

HOMO^LUMO (93%) 

c60ha 

941.0500 

0.3354 

HOMO^LUMO (88%) 

cSOha 

942.9108 

0.3294 

HOMO^LUMO (89%) 

c60he 

959.4006 

0.9232 

HOMO^LUMO (92%) 

c60hg 

959.8462 

0.6174 

HOMO^LUMO (93%) 

c54ha 

975.8625 

0.3409 

HOMO^LUMO (88%) 

c54hf 

985.0110 

0.9331 

HOMO^LUMO (92%) 

c74he 

994.1732 

1.0510 

HOMO^LUMO (87%) 

c58h32 

1002.9391 

0.3537 

HOMO^LUMO (87%) 

c74hd 

1010.7063 

1.0094 

HOMO^LUMO (82%) 

c74hc 

1024.2324 

0.3657 

HOMO^LUMO (85%) 

c62h34 

1026.2672 

0.3674 

HOMO^LUMO (85%) 

c66h36 

1045.7434 

0.3821 

HOMO^LUMO (84%) 

c74hb 

1053.2054 

0.3891 

HOMO^LUMO (82%) 

c70h38 

1062.0467 

0.3977 

HOMO^LUMO (82%) 

c60hc 

1064.8745 

1.0999 

HOMO^LUMO (92%) 

c74hg 

1191.0023 

1.1662 

HOMO^LUMO (92%) 

c60hd 

1193.2949 

0.9702 

HOMO^LUMO (91%) 

c74hf 

1209.8296 

1.2112 

HOMO^LUMO (91%) 

c74hh 

1359.3174 

1.0917 

HOMO^LUMO (91%) 

c74hi 

1405.2288 

0.9118 

HOMO^LUMO (89%) 

c74hw 

1793.4809 

1.0937 

HOMO^LUMO (90%) 

c74hz 

1794.7790 

1.1039 

HOMO^LUMO (90%) 
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c74hy 

1911.5531 

1.0413 

HOMO^LUMO (87%) 

c74hx 

1927.6016 

1.0427 

HOMO^LUMO (87%) 

c33hl5 

1980.5645 

0.0213 

HOMO(B)^L-bl(B) (61%), H-l(A)^LUMO(A) 
(13%) 

c74hp 

2049.6501 

0.9875 

HOMO^LUMO (84%) 

c74ho 

2068.8025 

0.9320 

HOMO^LUMO (84%) 

c74hn 

2104.2657 

0.8585 

HOMO^LUMO (83%) 

c74hm 

2176.2917 

0.7512 

HOMO^LUMO (83%), H-l^L-bl (10%) 

c78h24 

2427.7137 

0.2894 

HOMO^LUMO (14%), HOMO^L-tl (20%), 
H-l^LUMO (20%), H-l^L-tl (14%), H-l^L-b2 
(21%) 

c28hcc 

3158.0065 

0.2001 

HOMO^LUMO (87%) 

c46hl8 

4221.4279 

0.1481 

HOMO^LUMO (50%), H-l^L-bl (36%) 

c28hS 

4251.8291 

0.1646 

HOMO^LUMO (86%) 

c24hl4c 

4335.0817 

0.0951 

HOMO^LUMO (90%) 

c24hc 

4335.0817 

0.0951 

HOMO^LUMO (90%) 

c32he 

4359.4703 

0.1278 

HOMO^LUMO (90%) 

c28hM 

4413.7891 

0.1169 

HOMO^LUMO (88%) 

c32hc 

4659.2761 

0.1125 

HOMO^LUMO (90%) 

c28hR 

6156.0743 

0.0708 

HOMO^LUMO (89%) 

c74hj 

6515.1517 

0.0958 

HOMO^LUMO (58%), H-l^LUMO (34%) 

c28hss 

6687.3429 

0.0515 

HOMO^LUMO (91%) 

c38hk 

7505.0446 

0.0535 

HOMO^LUMO (89%) 

c28hj 

10818.7903 

0.0297 

HOMO^LUMO (92%) 

c30hj 

12167.1577 

0.0236 

HOMO^LUMO (92%) 

c34hg 

12251.3178 

0.0284 

HOMO^LUMO (92%) 

c26hR 

18126.2188 

0.0137 

HOMO^LUMO (93%) 


Table S2: System names, wavelength corresponding to, oscillator 
strength (Osc. Str.) of and molecular orbital contributions for 
^‘Pmax excitation” of all the GQDs whose p^ax is in IR-region. 

Names of the GQDs are as given in Ref 23 of the main article. 

All the corresponding structures can be obtained from this link. 
http://journals.aps.Org/prb/supplemental/10.1103/PhysRevB.81.085430/GNFs_PAHs_coord.tar.gz 


Table S3: System names, wavelength corresponding to “Pmas” ”Pmax’' 
excitations of all C28-GQDs whose p^ax is in IR-region are given. Values inside 
the paranthesis are the ZINDO/S results and the ones which are outside are the 
GAM-B3LYP/6-31-bg(d) results. 


System 

Pmax (nm) 

/3max (nm) 

c28hcc 

2966.1085 ( 3158.0065) 

436.1311 (460.6820) 

c28hj 

4125.9014 (10818.7903) 

346.5642 (429.1269) 

c28hM 

3035.8310 ( 4413.7891) 

434.3282 (460.0666) 

c28hR 

3145.1887 ( 6156.0743) 

393.7104 (454.2180) 

c28hS 

2896.1303 ( 4251.8291) 

425.8401 (491.0231) 

c28hss 

3408.0082 ( 6687.3429) 

397.7012 (448.3378) 
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Figure S2: Panel (a) shows the isosurfaces of HOMO of all C32 GQDs 
considered in this study. Panel (b) shows their corresponding absorption 
spectra. 


Table S4: System name, polarizability and 1st hyperpolarizability of 
C28-GQDs Values inside the paranthesis are the MOPAG results and the outside 
ones are CAM-B3LYP/6-31+g(d) results. 


System 

a (a.u) 

P (a.u) 

c28hcc 

744.480 (1502.36167) 

97762.095 (485305.1890) 

c28hM 

761.495 (1696.01550) 

199892.814 (1304022.4046) 

c28hss 

801.149 (1726.62126) 

437976.831 (2003517.2092) 

c28hj 

822.574 (2264.26242) 

869985.316 (5567249.5402) 
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Figure S3: Isosurfaces of HOMO of some of the C74 GQDs considered in this 
study. 
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Figure S4: Absorption spectra corresponding to the GQDs given in the above 
figure (Fig. S3). 
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